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Abstract 

We study some of the novel effects that arise when the QCD axion is placed 
in the "bulk" of large extra spacetime dimensions. First, we find that the mass 
of the axion can become independent of the energy scale associated with the 
breaking of the Peccei-Quinn symmetry. This implies that the mass of the ax- 
ion can be adjusted independently of its couplings to ordinary matter, thereby 
providing a new method of rendering the axion invisible. Second, we discuss the 
new phenomenon of laboratory axion oscillations (analogous to neutrino oscil- 
lations), and show that these oscillations cause laboratory axions to "decohere" 
extremely rapidly as a result of Kaluza-Klein mixing. This decoherence may 
also be a contributing factor to axion invisibility. Third, we discuss the role of 
Kaluza-Klein axions in axion-mediated processes and decays, and propose sev- 
eral experimental tests of the higher-dimensional nature of the axion. Finally, 
we show that under certain circumstances, the presence of an infinite tower of 
Kaluza-Klein axion modes can significantly accelerate the dissipation of the en- 
ergy associated with cosmological relic axion oscillations, thereby enabling the 
Peccei-Quinn symmetry-breaking scale to exceed the usual four-dimensional 
relic oscillation bounds. Together, these ideas therefore provide new ways of 
obtaining an "invisible" axion within the context of higher-dimensional theories 
with large-radius compactifications. 
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1 Introduction 



One of the pressing theoretical issues that confronts the Standard Model of par- 
ticle physics is to understand the origins of CP violation. While the weak inter- 
action is well-understood to lead to CP violation through complex phases in the 
Cabbibo-Kobayashi-Maskawa (CKM) fermion mass matrix, CP violation may also 
independently arise from the physics of the strong interaction. Ordinarily, one might 
have assumed the strong interaction to conserve CP. However, the U(l) problem [|TJ 
(namely, the inability to interpret the rj particle as the Nambu-Goldstone boson of a 
spontaneously broken U(l) flavor symmetry) turns out to require a non-trivial QCD 
vacuum structure which in turn naturally leads to CP violation. Specifically, 
one finds that the effective Lagrangian describing the strong interaction should be 
augmented by an additional CP- violating contribution 

£ eS = £qcd + e J^FfF^ (1.1) 

where the parameter receives two contributions: 

= + ArgdetM . (1.2) 

Here is the strong-interaction 0-angle reflecting the non-trivial nature of the QCD 
vacuum, and Argdet M (with M denoting the CKM matrix) is the contribution 
arising from weak interactions. However, measurements of the neutron electric dipole 
moment place the stringent experimental bound 

< 10~ 9 . (1.3) 

Explaining the small size of is the strong CP problem pi. 

The most elegant explanation of the strong CP problem is the Peccei-Quinn (PQ) 
mechanism ||, i n which is set to zero dynamically as the result of a global, sponta- 
neously broken £7(1) Peccei-Quinn symmetry. However, associated with this symme- 
try is a new Nambu-Goldstone boson, the axion |5|, |6| , which essentially replaces the 
parameter in the effective Lagrangian. This then results in an effective Lagrangian 
of the form || 

£ eff = £qcd + \d,ad»a + ^JPfFF^ (1.4) 

Jpq o2n z 

where /pq is the axion decay constant, associated with the scale of PQ symmetry 
breaking. Here £ is a model-dependent parameter describing the PQ transformation 
properties of the ordinary fermions, and we have not exhibited further Lagrangian 
terms which describe axion/fermion couplings. The mass of the axion is then expected 
to be of the order 

A 2 

m a ~ -f — (1.5) 
JPQ 
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where Aqcd ~ 250 MeV; likewise, the couplings of axions to fermions are suppressed 
by a factor of 1//pq. Thus, heavier scales for PQ symmetry breaking generally imply 
lighter axions which couple more weakly to ordinary matter. 

Ordinarily, one might have liked to associate the scale /pq with the scale of 
electroweak symmetry breaking, implying an axion mass near the electron mass 
m a « O(10 2 ) keV. However, experimental searches for such axions have so far been 
unsuccessful 0, and indeed only a narrow allowed window currently exists 

10 8 GeV < /p Q < 10 12 GeV. (1.6) 

This then implies an axion which is exceedingly light 

l(T 5 eV < m a < lO^eV , (1.7) 

and whose couplings to ordinary matter are exceedingly suppressed. These bounds 
generally result from various combinations of laboratory, astrophysical, and cosmo- 
logical constraints. In all cases, however, the crucial ingredient is the correlation 
between the mass of the axion and the strength of its couplings to matter, since both 
are essentially determined by the single parameter /pq. 

In this paper, we point out that this situation may be drastically altered in the- 
ories with large extra spacetime dimensions. Since their original proposal P], such 
theories have recently received considerable attention because of their prospects for 
lowering the fundamental GUT scale 0, the fundamental Planck scale [110], and the 



fundamental string scale As a result of these developments, it is now under- 

stood that all three of these scales may be adjusted to arbitrary values, perhaps even 
values in the TeV-range, without violating experimental constraints. Of course, if 
one lowers these fundamental scales below the PQ symmetry- breaking scale in (|1.6| ). 
then it might seem difficult to preserve the axion solution to the strong CP problem. 



In fact, this observation has been used [12| to argue that the fundamental scales of 



physics should be taken at some intermediate scale near the PQ scale. However, this 
argument assumes that the PQ mechanism itself remains untouched by the presence 
of the extra large dimensions. 

In this paper, we shall generalize the PQ mechanism to higher dimensions. More 
specifically, we shall consider the consequences of placing the PQ axion in the "bulk" 
(i.e., perpendicular to the p-brane that contains the Standard Model) so that the 
axion accrues an infinite tower of Kaluza-Klein excitations. Placing the QCD axion 



in the bulk has also been discussed previously in Refs. |L0L O], and is completely 



analogous to recent ideas concerning the placement of the right-handed neutrino in 



the bulk PH, |I5| in order to lower the neutrino seesaw scale and obtain neutrino 
oscillations with light (or even vanishing [fPfl ) neutrino masses. 

As we shall see, placing the axion in the bulk has a number of surprising con- 
sequences for axion physics, all of which provide new ways of rendering the axion 
"invisible" . 
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• First, in Sect. 2, we shall show that the mass of the axion can be decoupled from 
the PQ symmetry-breaking scale, and under certain circumstances is essentially 
set by the radius R of the bulk: 

m a < 0(R- r ) . (1.8) 

Thus, remarkably, the allowed range ( |1.7|) for the axion mass is consistent with 
radii in the millimeter or sub-millimeter range! More importantly, however, 
this result implies that it is possible to adjust the PQ symmetry breaking scale 
independently of the axion mass in order to control the strength of the couplings 
of this axion to ordinary matter. This might be useful in order to provide an 
alternative explanation for an "invisible" axion. 

• Second, in Sect. 3, we shall show that the presence of an infinite tower of 
Kaluza- Klein axion modes can induce the novel phenomenon of laboratory axion 
oscillations . These oscillations are completely analogous to laboratory neutrino 
oscillations, but we shall find that under certain circumstances they lead to a 
complete and rapid decoherence of the axion field. This implies that an axion, 
once produced in the laboratory, will "decohere" extremely rapidly. This is 
therefore a second higher-dimensional phenomenon that can contribute to the 
"invisibility" of the axion under certain circumstances. 

• Third, in Sect. 3, we shall also discuss the role played by the excited axion 
Kaluza-Klein states in axion-mediated processes and decays, and propose sev- 
eral experimental methods of detecting their existence. These would therefore 
provide direct experimental tests of the higher- dimensional nature of the axion 
field. 

• Finally, in Sect. 4, we shall show that under certain circumstances, the pres- 
ence of the infinite tower of Kaluza-Klein axion modes can actually accelerate 
the dissipation of the energy associated with cosmological relic axion oscilla- 
tions. This implies a weakening of the usual four- dimensional cosmological 
relic oscillation bounds on the PQ symmetry-breaking scale, which in turn may 
permit axion/matter couplings to be suppressed even more strongly than in 
the usual four- dimensional case. This is therefore a third higher-dimensional 
feature leading to an "invisible" axion. 

Together, these results suggest that "invisible" axions can emerge quite naturally 
within the context of higher-dimensional theories with large-radius compactifications, 
and have significantly different phenomenologies than they do in four dimensions. 
Moreover, many of our results apply to bulk fields in general, and transcend the 
specific case of the axion. This is illustrated in Sect. 6, where we consider the conse- 
quences of extra dimensions for another hypothetical particle, the so-called Standard- 
Model dilaton. 
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2 A higher-dimensional Peccei-Quinn mechanism 



In order to generalize the Peccei-Quinn (PQ) mechanism to higher dimensions, 
we will assume that there exists a complex scalar field in higher dimensions which 
transforms under a global £7(1 )pq symmetry: 

-> e iA <f) . (2.1) 

This symmetry is assumed to be spontaneously broken by the bulk dynamics so that 
(0) — /pq/ \/2) where fpq is the energy scale associated with the breaking of the PQ 
symmetry. We thus write our complex scalar field in the form 

/PQ Ja/fpQ 



V2 



e «v/pq (2.2) 



where a is the Nambu-Goldstone boson (axion) field. If we concentrate on the case 
of five dimensions for concreteness, then the kinetic-energy term for the scalar field 
takes the form 

5 k .e. = J d 4 xdy M s d M <j)*d M <f) = J d^xdy M s \d M ad M a (2.3) 

where M s is a fundamental mass scale (e.g., a Type I string scale), and where we have 
neglected the contributions from the radial mode. Here x M are the coordinates of the 
four uncompactified spacetime dimensions, y is the coordinate of the fifth dimension, 
and the M spacetime index runs over all five dimensions: x M = (x^,y). Note that 
there is no mass term for the axion, as this would not be invariant under the {7(1)pq 
transformation 

a - a + / PQ A . (2.4) 

Furthermore, as a result of the chiral anomaly, we will also assume a bulk/boundary 
coupling of the form 



'coupling 



J d ± xd y-Lj?_ a Fr p^ a 5(y) (2.5) 



fpq 32tt2 



where F^ ua is the (four-dimensional) QCD field strength describing the QCD gauge 
fields which are confined to a four- dimensional subspace (e.g., a D-brane) located at 
y — 0, and where £ is a model-dependent quantity parametrizing the strength of the 
axion couplings to matter. Thus, our total effective five-dimensional axion action 
takes the form 



?eff 



= J d A x dy 



\M s d M ad M a + A_^aFrP^ a 5(y) 



(2.6) 



While we have assumed that the spontaneously broken C/(1)pq is parametrized 
by fpQ, one still has to address the fact that gravitational effects can also break the 
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L r (l)pQ symmetry. In other words, gravitational interactions do not respect global 
symmetries. Ultimately, this can lead to a gravitational contribution to the axion 
mass which is not necessarily suppressed. In this paper, however, we will assume 
that the gravitational contributions to the axion mass are indeed suppressed, and 
that U(1)pq remains a valid symmetry even in the presence of gravitational effects. 
For example, such suppression might arise due to the suppression of gravitational 
interactions across a large bulk, discrete symmetries of the sort that might come 
from Scherk-Schwarz compactifications, or other large-radius effects. 

In order to obtain an effective four-dimensional theory, our next step is to com- 
pactify the fifth dimension. For simplicity, we shall assume that this dimension is 
compactified on a Z 2 orbifold of radius R where the orbifold action is identified as 
y — ► —y. This implies that the axion field will have a Kaluza-Klein decomposition of 
the form 



a[x 



y) = £ a n (x») cos ( ^ ) (2.7) 



R 



n=0 

where a n (x^) G 1R are the Kaluza-Klein modes and where we have demanded that 
the axion field be symmetric under the Z 2 action (in order to have a light zero-mode 
that we can identify with the usual four-dimensional axion). 

In principle, we should also allow for the possibility that the axion field winds 
non-trivially around the extra compactified dimension, with winding number w. This 
possibility of winding arises because a is really only an angular variable, as evident 



from (2.2), and would imply that we should introduce an extra term wf^qy/R into 



the mode expansion ( |2.7[ ). Note, in particular, that such windings cannot be removed 



by global Peccei-Quinn transformations of the form ( |2.4| ). However, such a winding 
term would not be invariant under the orbifold symmetry y — > —y, and thus only 
the unwound configuration w = survives the orbifold projection. Even if we were 
to compactify the axion field on a circle rather than an orbifold, such a term would 
contribute only an overall additive constant to the resulting effective potential for 
the axion field, and would not change any of the subsequent physics. It is therefore 
sufficient to restrict our attention to the unwound configuration with w — 0. 

It is also interesting to note that with respect to the four-dimensional Kaluza- 
Klein axion modes a n , the Peccei-Quinn transformation (|2.4j) takes the form 



«o — > a o + /pqA ^2 g\ 

a-k — ► a-k for all k > . 

Thus, we see that only ao serves as the true axion transforming under the PQ trans- 
formation, while the excited Kaluza-Klein modes remain invariant. 

Substituting (|2.7|) into (|2.6|) and integrating over the fifth dimension, we then 
obtain an effective four-dimensional Lagrangian density 

oo oo 2 f a 2 / 00 \ 

C cS = £qcd + |E(W - I E m«n + T- 4^ E r " a " (2.9) 

n=0 n=l /PQ l)Z7r \n=0 / 
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where 

f 1 if n = 



r » 3 1^2 if „ > . < 210 » 



Note that in order to obtain Q2.9| ), we must individually rescale each of the Kaluza- 
Klein modes a n in order to ensure that they have canonically normalized kinetic- 
energy terms. It is this that produces the relative rescaling coefficients r n in fl2.10|) . 
We have also defined /pq = (VM s ) 1 / 2 / P q, where V is the volume of our compactified 
space. For 8 extra dimensions, this definition generalizes to /pq = (VM$ ) 1 ^ 2 /pq. For 
example, assuming a 5- dimensional toroidal compactification implies V = (2ttB) s , 
resulting in the relation 

/pq = (2nRM s ) s / 2 / PQ . (2.11) 

Note that while /pq sets the overall mass scale for the breaking of the Peccei- 
Quinn symmetry, it is the volume-renormalized quantity /pq that parametrizes the 
coupling between the axion and the gluons. In general, since M s Br 1 , we find that 
/pq ^> /pq- Therefore, as pointed out in Ref. [[H]], this volume-renormalization of 
the brane/bulk coupling can be used to obtain sufficiently suppressed axion/gauge- 
field couplings even if /pq itself is taken to be relatively small. In other words, even 
if we demand that /pq be in the approximate range 10 8 GeV < /pq < 10 12 GeV, 
the fundamental Peccei-Quinn symmetry-breaking scale /pq can be substantially 
reduced, potentially all the way to the TeV-range. This volume suppression is thus 



one higher-dimensional way |IU| of avoiding the need for a high fundamental Peccei- 
Quinn scale /pq. 

Two comments are important at this stage. First, in toroidally compacti- 
fied higher-dimensional theories with reduced Planck scales, the higher-dimensional 
Planck scale M* and the usual four-dimensional Planck scale Mpi anc k are related to 
each other via [|10j 

M Planck = (27T-RM*)"/ 2 M* (2.12) 

where n is the total number of extra spacetime dimensions in the bulk. Identifying 
M* ~ M s , we therefore see that the n-dimensional volume factor (2irBM s ) n / 2 must 
already be adjusted in order to account for the difference between M s ~ C(TeV) 
and M Planck ~ (9(10 19 GeV). If we were to take 6 = n for the current axion case, 
this would imply either that /pq ~ Mpi anc k (which would presumably overdose the 
universe), or that M s -C O(TeV) (which would clearly violate current experimental 
bounds). Therefore, if we assume an isotropic compactification with all radii taken 
equal, we see that an intermediate scale /pq can be generated via ( [2.1 1| ) only if we 
have 5 < n. In other words, under these assumptions, the axion must be restricted 
to a subspace of the full higher-dimensional bulk. This has already been pointed 
out in Ref. ||13|| , and is analogous to similar restrictions that arise in the case of 



higher-dimensional neutrinos [14, 15 



Our second comment concerns the relevance of the mass scale /pq that is generated 
by this volume factor. Of course, it is apparent from ( p.9|) that /pq (rather than /pq) 
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sets the scale for couplings between gauge fields and individual axion modes a^. 



However, we have also seen in ( |2.9| ) that the gauge fields couple not to an individual 



axion mode a^, but rather to the linear superposition 

^max | / ^max \ 

a' = —= r na n = -7= a + v / 2^ a « ( 2 - 13 ) 
V-/V n=0 V i\ \ n= i ) 

where 

N = l + 2n max . (2.14) 

Here n max is a cutoff, determined according to the underlying mass scale M s (which 
sets the limit of validity of our higher-dimensional effective field theoryQ). Taking 
n max « RM S ^> 1, we then find that the axion/gluon coupling in ( |2.9| ) takes the form 



N , ~ 1 , 

f a Fa Ffj.ua ~ j— f a, F£ F^a . (2-15) 

/PQ V^JPQ 

Thus it is actually /pq, rather than /pq, that sets the scale for axion couplings 
involving the entire Kaluza-Klein linear superposition a'. In other words, the effects 
of the volume factor in (|2.11 ) are cancelled by the normalization of the Kaluza-Klein 



linear superposition a'. Of course, this is expected from the perspective of the higher- 
dimensional theory in which /pq is the only fundamental mass scale, and indeed 
this cancellation of the volume factor persists for any number of extra spacetime 
dimensions. Because we expect /pq /pq in scenarios with large extra spacetime 
dimensions, we see that axion couplings involving the linear superposition a' are 
relatively strong, and pose a serious threat to the invisibility of the higher-dimensional 
axion. We shall discuss how this problem may be overcome in subsequent sections. 

Similar observations also apply for axion couplings to Standard-Model fermions. 
The invariance under the C/(1)pq transformation a — > a + /pqA implies that axions 
can be at most derivatively coupled to fermions carrying a PQ charge. If we assume 
that these fermions are also restricted to the D-brane at y = (as would be the case 
for all Standard- Model fermions), then this axion/fermion coupling is restricted to 
take the form 

S a w ~ f d 4 x J-id^o)^^) . (2.16) 



Here a\ y=0 is the full five-dimensional bulk axion field evaluated at y = 0, and the 
bulk/brane coupling strength /pq is defined in ( |2.11|) . Note that /pq is the same 



volume-rescaled axion decay constant that parametrizes the couplings to the gauge 
fields, with the volume factor emerging just as for the axion/gauge couplings. Using 
the Kaluza-Klein decomposition (|2.7| ), the action then becomes 



DC 



Saw ~ — E r ™ /rf 4 x(^a n )(^7 5 ^) ~ — /^(VHVWV) > ( 2 - 17 ) 

n=0 J 



* Note that the heavy modes with masses of order M s can be treated within field theory along 
the lines discussed in Ref. jl6|], and within string theory along the lines discussed in Ref. p7f. 
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and we see that once again the entire Kaluza-Klein linear superposition a' couples to 
the charged fermions. This is completely analogous to the situation in ( |2.9|) for the 
gauge fields. It is important to stress that this need not have been the result from a 
purely four- dimensional perspective. Indeed, given the Peccei-Quinn transformation 
properties (|2.8j), we see that it is only the zero-mode ao which requires a derivative 
coupling to fermions; the other axion modes are a priori free to have non-derivative 
couplings. It is therefore only the higher- dimensional structure of the axion field 
that forces all of the axion modes to have identical derivative couplings to charged 
fermions. Moreover, we see from (p. 17 ) that while the mass scale for the couplings 



of individual Kaluza-Klein axion fermions to fermions is set by /pq, the mass scale 
for the coupling of the full linear superposition a! to fermions is set by /pq. 

Let us now proceed to verify that this higher-dimensional Peccei-Quinn mecha- 
nism still cancels the CP-violating phase, and use this to calculate the mass of the 
axion. Given the Lagrangian (|2.9|) , we observe that a' ~ X^o r " a « serves as the 
overall quantity that parametrizes the size of CP symmetry breaking. Applying the 
one-instanton dilute-gas approximation, it is straightforward to show that 

(Fr^ a ) = - A QCD sin (J- f) r n a n + &] . (2.18) 

This gives rise to an effective potential for the axion modes in the QCD vacuum 

00 n 2 a 2 ( £ 00 -Y 

In order to exhibit the Peccei-Quinn mechanism, we now minimize the axion effective 
potential, 

dV n 2 

da--& an + rn S^ QCDDlli V^o 
yielding the unique solution 

(a ) = ^(-6 + fvr) , £ G 22 

(a*) = for all k > . (2.21) 

Note that while any value i G 2 provides an extremum of the potential, only the 
values i G 2Z provide the desired minimum of the potential. Thus, this higher- 
dimensional Peccei-Quinn mechanism continues to solve the strong CP problem: we 
see that ao is the usual Peccei-Quinn axion which solves the strong CP problem by 
itself by cancelling the angle, while all of the excited Kaluza-Klein axions for 
k > have vanishing VEVs. This makes sense, since only ao is a true massless 
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Nambu-Goldstone field from the four-dimensional perspective of ( |2.9| ). This is also 
evident from the Peccei-Quinn transformation properties ( |2.8|) . 

However, as we shall now show, these excited Kaluza-Klein axion states never- 
theless have a drastic effect on the axion mass matrix. In order to derive the mass 
matrix, we consider the local curvature of the effective axion potential around its 
minimum: 

d 2 v 



From this we obtain 

2 



M 2 , = 

nn da, n da n/ 



n 2 A 4 
C 2 9 A QCD 
+ f 777^7 — r nTn< COS 



(a) 



327T 2 



./l 2 



PQ 



JpQ 



r n a n + 9 



n=0 



(2.22) 



(2.23) 



and in the vicinity of the minimum ( 2.21|) this becomes 

.2 9 2 



Mt 



n 2 

W 6n 



+ e 



,2 A 4 



QCD 



32n 2 f 2 



PQ 



Let us now define 



m 



PQ 



y = 



n 2 A 4 
-2 9 iv QCD 



32tt 2 



PQ 



m PQ R 



(2.24) 



(2.25) 



Thus mpQ is the expected mass that the axion would ordinarily have taken in four 
dimensions (depending on /pq rather than /pq itself), and y is the ratio of the scale 
of the extra dimension to mpq. Our mass matrix then takes the form 



or equivalently 



M' 



M' 



m 



PQ 



— m PQ {j"n r n' 


i 2 2 

+ y n 






1 V2 


V2 


y/2 




V2 2 + y 2 


2 


2 




V2 2 


2 + % 2 


2 




y/2 2 


2 


2 + 9y 2 





\ 



(2.26) 



^2.27) 



Note that the usual Peccei-Quinn case corresponds to the upper-left lxl matrix, 
leading to the expected result .M 2 = mpq. Thus, the additional rows and columns 
reflect the extra Kaluza-Klein states, and their physical effect is to pull the lowest 
eigenvalue of this matrix away from mpq. 
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Deriving the condition for the eigenvalues of this matrix is straightforward. Let 
us denote the eigenvalues of this matrix as A 2 rather than A because this is a (mass) 2 
matrix. We then find that the eigenvalues are given as the solutions to the transcen- 
dental equation 

7rA / 7rA\ ~ 9 

— cot — = A 2 2.28 

v V v J 

where we have defined the dimensionless eigenvalue 

A = A/mp Q . (2.29) 

In terms of dimensionful quantities, this transcendental equation takes the equivalent 
formp] 

7ri?A cot( 7 r J RA) = -4~ . (2.30) 

For each eigenvalue A, the corresponding normalized mass eigenstate a\ is exactly 
given by 

oo 

a\ = u ^kCik (2.31) 

fc=0 



where are the Kaluza-Klein axion modes given in (|2.7| ) and where U\k is the unitary 
matrix that diagonalizes M 2 . This matrix is given by 

u » s {>rh) Ax (2 - 32) 

-4a = ^ (X' + l + ^/y 1 )' 1 ' 2 . (2.33) 

Note that the unitarity of the matrix U implies that J2x I^ao| 2 = 1? which in turn 
implies 

E A l = 1 ■ ( 2 - 34 ) 



where 



t Interestingly, this eigenvalue equation is identical to that which emerges when the right- 
handed neutrino vr is placed in the bulk, with the mass scale mpQ in the axion case corresponding 
to the Dirac coupling m in the neutrino case. This implies that there is a formal relation between the 
Kaluza-Klein axion modes and the Kaluza-Klein neutrino modes. Remarkably, this correspondence 
exists even though the axion and right-handed neutrino have different spins, and even though the 
mechanisms for mass generation are completely different in the two cases. Moreover, in Sect. 6, 
we shall demonstrate that the same mass matrix and eigenvalue equation also emerge when the 
Standard-Model dilaton is placed in the bulk of extra spacetime dimensions. This suggests that 
many of the higher-dimensional phenomena to be discussed in this paper (such as laboratory and 
cosmological relic axion oscillations) may have a correspondingly general phenomenology that is 
equally applicable to neutrinos, dilatons, as well as other bulk fields that transform as singlets 
under the Standard-Model gauge group. 



10 



For future reference, we also record another useful identity which will be proven in 
Sect. 4: 

£AM* = 1. (2-35) 

A 

Finally, combining ( 2.28Q and Q2.32 ), it is straightforward to show that 



E r kUxk = A 2 A\ . (2.36) 



fc=0 



Note that all three of these identities hold for all values of y = (mpqR) 1 . This in 
turn allows us to rewrite fl2.13p in the form 



1 E A 2 ^a«a • (2.37) 



N A 

We can check that ( p.30| ) makes sense in the limit R — > 0. In this limit, the 
Kaluza-Klein states become infinitely heavy and decouple; thus we should be left 
with the lightest eigenvalue A = mpq. And this is indeed what happens: as R — > 0, 
we see tiRX cot(nR\) — > 1, whereupon we obtain A = mpq, with all other eigenvalues 
infinitely heavy. As R becomes larger, the effect of the extra large dimension is felt 
through a reduction of this lowest eigenvalue. Thus, as R increases, the mass of the 
lightest axion decreases. 

One important consequence of ( |2.30| ) is that the lightest axion mass eigenvalue 
m a is strictly bounded by the radius 



m n 



< fiT 1 . (2.38) 



This result ( |2.38| ) holds regardless of the value of mpq. Thus, in higher dimensions, 
we see that when mpq > the size of the axion mass is set by the radius R and 

not by the Peccei-Quinn scale fpq. In Fig. [I] we show the value of the axion mass 
m a as a function of y 1 = m^qR. Of course, for mpQR —>■ (corresponding to either 
R — >• or mpQ — > 0), we see that we indeed have the expected result m a ~ mpq. 
This is indicated by the diagonal dashed line in Fig. 0. However, as mpqR increases, 
we see that the axion mass departs from this expected linear behavior, and instead is 
bounded by the inverse radius of the extra spacetime dimension. In fact, from Fig. |], 
we see that we can approximate the mass of the axion as 



min (fiT 1 , m PQ ) . (2.39) 



Thus the mass of the axion is determined solely by the radius of the extra spacetime 
dimension when f-R -1 ^ ^pq- 

As a result of this unexpected higher- dimensional behavior for the axion mass, we 
see that when mpq > \R~ l i the Peccei-Quinn scale fpq essentially decouples from 
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Figure 1: The mass of the axion zero-mode as a function of the dimensionless product 
mpQ.R = where w-pq ~ Aqcd//pQ- Although this reproduces the expected result 
m a mpQ when mpq <C R^ 1 (corresponding to the diagonal dashed line \R = mpQR), 
we see that the axion mass is strictly bounded by the inverse radius in higher dimensions, 
with the precise value of mpQ essentially decoupling for topq > This implies that 

in higher dimensions, the Peccei-Quinn scale /pn can be adjusted arbitrarily within this 
parameter range without affecting the axion mass. 



the axion mass! Indeed, as long as m P q > \R~ 1 , we see that m a < \R~ l regardless 
of the specific sizes of mpQ or Aqcd- 

This observation has a number of interesting implications. First, given (|2.38| ), we 
see that an axion mass in the allowed range ( |1.7| ) is already achieved for R in the sub- 
millimeter range, independently of mpq! This therefore provides further motivation 
for such submillimeter extra dimensions. 

Second, and even more importantly, we now have the surprising result that mpQ 
can still be lowered or raised arbitrarily without upsetting the constraint ( |1.7| ), pro- 
vided m P q > \R~ X - In other words, having already satisfied the axion mass con- 
straints by appropriately choosing the value of R, we are now essentially free to tune 
mpQ (or equivalently the fundamental Peccei-Quinn symmetry breaking scale /pq) 
in such a way as to weaken the axion couplings to matter to whatever values are 
required to make the axion sufficiently invisible. This may therefore provide a new 
method of obtaining an invisible axion. 
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3 Laboratory axion oscillations 



In this section we discuss the novel possibility of laboratory axion oscillations. 
Note that these are not the traditional cosmological relic axion oscillations (which 
will be discussed in Sect. 4), but rather "laboratory" axion oscillations that now arise 
in higher dimensions because the physical Peccei-Quinn axion ao is no longer a mass 
eigenstate. These laboratory axion oscillations are therefore completely analogous 
to laboratory neutrino oscillations, which similarly arise because the neutrino gauge 
eigenstates differ from the neutrino mass eigenstates. In the present axion case, such 
oscillations are possible because of the non-diagonal nature of the axion mass matrix 
( |2.27D . Therefore, as it propagates, it is possible for the four- dimensional axion zero- 
mode ao to oscillate into any of the higher-frequency axion modes in its Kaluza-Klein 
tower. Indeed, as we shall see, these oscillations can provide yet another mechanism 
that may contribute to the "invisibility" of the axion. 

3.1 Laboratory oscillations of the axion zero-mode ao 

It is straightforward to calculate these oscillation probabilities in terms of the 
mass mixing matrix ( 2.27 ) and the [/-matrix ( |2.32j) that diagonalizes it. From a 



four-dimensional perspective, we see from fl2.8|) that only the zero-mode ao serves as 
a bona-fide axion transforming under the PQ transformation. Therefore, let us first 
calculate the probability that this four- dimensional axion a oscillates into any of its 
corresponding Kaluza-Klein excitations as it propagates, or conversely the probability 
that the four- dimensional axion ao is preserved as a function of time. Assuming 
that the axion is given an initial highly relativistic momentum p, we find that the 
probability that ao is preserved is given by 



]T^e- r ^ 2 exp 



2p 



X^e" 1 * + 2£ A^e-^+ r v)«/a cos ( \^L . (3.1) 



A'<A 



2p 



Here A\ is defined in (|2.33| ), and T\ is the decay width of the corresponding Kaluza- 
Klein axion. Thus, even though Pq^o = 1 at the initial time t = 0, we see that this 
probability decreases at later times and ultimately oscillates around a period-averaged 
value 

(*W*)> = E^- r ^ (3-2) 

A 

which itself diminishes exponentially with time. Note that the calculations leading 
to these results are similar to the higher-dimensional neutrino oscillation calculations 
in Ref. |4|. 



It is important at this stage to separate two effects which influence the axion 
preservation probability. The first is the oscillation itself, which arises due to the 
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non-trivial axion mass matrix and which reflects the mixing of the excited Kaluza- 
Klein axion states. It is this oscillation which is our focus in this section. By contrast, 
the second effect is axion decay as reflected in the decay widths I\. In general, these 
decay widths result from the dominant decay mode d\ — > 77, and therefore scale 
as T\ « A 3 //|q. This implies that for the lowest Kaluza-Klein eigenvalues and 
sufficiently small times, the product F\t is typically extremely small. This is in 
accordance with our expectation that the usual four-dimensional axion is extremely 
stable. For example, taking Ao ~ mpQ ~ 10~ 5 eV, we see that T t <C 1 for all times 
t < 10 48 seconds. This upper limit exceeds the age of the universe by 20 orders of 
magnitude. Therefore, particularly for the lowest Kaluza-Klein eigenvalues, it is safe 
to neglect these decay widths entirely, and concentrate solely on the oscillations. For 
example, since A\ decreases rapidly as a function of A, only the lowest eigenvalues 
dominate the sum in (|3.2|) . We therefore find that we can approximate 

(W*)> « E4. (3.3) 

A 

which is independent of time. We stress, however, that the decay widths T\ grow 
rapidly as a function of the mass of the Kaluza-Klein eigenstate a\. We have therefore 
included these decay widths in ( |3.1| ) and ( pT2| ) for completeness, and will discuss the 
effects that they induce more carefully at the end of this section. 

In Fig. 0(a), we have plotted the behavior of P ^ (t) as a function of time, taking 
a reference value y = (mpQ_R) _1 = 0.4. Note that the jaggedness of the probability 
curve reflects the multi-component nature of the oscillation in which many different 
individual Kaluza-Klein oscillations interfere with incommensurate phases. Although 
this oscillation clearly leads to both axion deficits and axion regenerations, we see 
that while the axion regenerations are nearly total, the axion deficits are not total. 
This is in marked constrast to the results from a simple two-state oscillation. We 
also observe that these oscillations are approximately periodic, with a wavelength 
set by the lowest-lying eigenvalue difference. This is because it is the lowest-lying 
Kaluza-Klein axions that play the dominant role in producing this oscillation. As a 
result of this fact, we see that (Po^o(^)) is effectively constant as a function of time, 
in accordance with ( |3.3| ). Indeed, the interpretation of this oscillation is completely 
analogous to that given in Ref. |14| for higher-dimensional neutrino oscillations. We 



shall discuss the possibilities for experimentally detecting such oscillations at the end 
of this section. 

In Fig. §(b), we have plotted the time-averaged probability (P ^ (t)) given in 
(|3.3|) as a function of y = (mpQ-R) -1 . For y ^> 1 (corresponding to the usual four- 
dimensional limit), we see that (Po^o(^)) — > 1, as expected, reflecting the fact that our 
single axion field cannot oscillate because its Kaluza-Klein states are infinitely heavy 
and essentially decouple. More interestingly, however, we see that in the y <C 1 limit 
(corresponding to extremely large radii or equivalently a quasi-continuous spectrum 
of light Kaluza-Klein modes), the oscillation probability ( |3.3| ) approaches a fixed 
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20 40 60 1 2 3 4 

(m a PQ /2p) t y = (m^R)- 1 

Figure 2: Higher-dimensional axion oscillations, as discussed in the text, (a) The probability 
Po^o(t) that the axion zero- mode is preserved as a function of time, taking a reference 
value y = (mpQ.R) _1 = 0.4. (b) The period-averaged preservation probability (Pq-»o(*)) 
as a function of y. The limit )/ > 1 corresponds to the usual four-dimensional case, while 
the opposite y <C 1 limit corresponds to an extremely large extra dimension with very light 
Kaluza-Klein states. 



value 

lim (P _o(t)) = lim E^A = \ ■ (3.4) 

Thus, for extremely large radii, we expect to see on average only 2/3 of the axion 
flux that would have appeared in the four- dimensional case. 

3.2 Laboratory oscillations of the axion superposition a' 

For many practical purposes, there exists a different probability that may be 
more relevant as a measure of laboratory axion oscillations. As we have seen in 
Sect. 2, Standard-Model gauge bosons and fermions generically couple not to a , 
but rather to the linear superposition a' given in ( |2.13| ). Thus, in any laboratory 
process that produces axions or is mediated by axions, a more crucial oscillation 
probability is the probability P a >^ a >(t) that this particular linear combination a' is 
preserved as a function of time. Indeed, we have already seen in Sect. 2 that while 
the Standard- Model couplings to individual axion modes scale as 1//pq and hence 
are already somewhat "invisible", the couplings to a' scale as 1//pq and hence are 
significantly larger. Such couplings therefore pose the largest immediate threat to 
axion invisibility. 
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The calculation of the probability P a /_> /(t) proceeds in an analogous manner. In 
general, the amplitude for an axion transition a k — > ae is given by 



A k ^{t) = J2U xe U* k e 



-V x t/2 „-i\ 2 t/2p 



(3.5) 



where U\k are the (real) unitary matrix elements defined in ( |2.32| ) and where we have 
omitted an overall (k, £)-independent phase. The probability that a' is preserved as 
a function of time is therefore given by 

2 



Pa'^a 



r k r e A k _i(t) 



k,£=0 



(3.6) 



where the normalization factor N is given in (gig ). In evaluating (gj|) , it is con- 
venient to use the identity (|2.36| ) in order to perform the Kaluza-Klein summations. 
We thus find 



1 

AT2 



2 £ A 4 (A') 4 ^^e-^ +r V)*/ 2 cos 

A'<A 



[A 2 -(AQ 2 ]r 
2p t 



implying a period-averaged probability 

(Pa>-+a>(f)) 



N 2 



(3.7) 



(3. 



Just as with Po^o(^), we will find it convenient to distinguish between two different 
effects: the overall "damping" that arises due to axion decays (encoded within the 
decay widths I\), and the oscillations that arise due to the non-trivial mixings of 
the excited Kaluza-Klein states. In order to concentrate on the latter effect, we shall 
therefore set T\ = for simplicity. Moreover, as we shall see, this assumption will 
not change our phenomenological results regardless of the time interval t in (|3.10|) . 
Thus, taking T\ = 0, we see that ( |3.7| ) reduces to 



EA 8 4 + 2 E A^TAlAl cos 

A A'<A 



[A 2 - (A') 2 ]t 
2p 



implying a time-averaged probability 

{Pa^a'it)) 



N 2 



(3.9) 



(3.10) 



It is straightforward to evaluate this time-averaged probability for different values 
of y. Remarkably, however, we find that 



lim (P^(t)) = ! 



(3.11) 
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We therefore have virtually no probability for detecting the linear combination a' 
at any later time after it is produced! Indeed, at the initial time t = 0, the axion 
probability starts at 1. This is guaranteed by the unitarity of the [/-matrix in ( |2.32| ), 
and may be verified directly from (|3.9| ). However, for t > 0, the multi-component 
Kaluza-Klein oscillations drive the net probability rapidly to zero. Indeed, as we 
shall see, this holds for all i > 0, and does not rely on taking the t — > oo limit. At no 
later time does a macroscopic axion regeneration appear. It is for this reason that it 
is justified to set T x = in ( |3.7| ) and ( P-8|) . 




0.001 0.002 0.003 0.004 
(nip /2p) t 



Figure 3: The axion preservation probability P a >^a'(t) 

as a function of the number 7i max 
of Kaluza-Klein states which are included in the system. For this plot we have set y = 15, 
and taken (a) n max = 1; (b) n max = 2; (c) n max = 3; (d) n max = 5; and (e) n max = 30. 
As n max increases, the axion probability rapidly falls to zero as a result of the destructive 
interference of the Kaluza-Klein states, and remains suppressed without significant axion 
regeneration at any later times. Note, in particular, that the "spikes" in this plot are also 
suppressed as ~ l/n max , and vanish for large n max . Thus, for significantly large n max , the 
destructive interference of the excited Kaluza-Klein states causes the axion to "decohere" , 
implying that there is negligible probability for subsequently detecting the original axion 
state at any future time. 

Of course, these conclusions rely on taking the n max — > oo limit in ( |3.11| ). However, 
even if we truncate n max at a finite value ~ 0(M string R) (reflecting the expected limit 
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of validity for our effective field-theoretic treatment), we find 

(iV+«'(*)> S (M^gi?)' 1 . (3.12) 

For M string « 10 TeV and R ~ 1 millimeter, this implies 

(P.w(*)> S 0(KT 17 ) • (3.13) 

Such an axion-mediated process therefore continues to be truly invisible in the sense 
that the a' state literally "disappears" over a very short time interval, with only a 
vanishingly small probability for detection of the a' state at any later time. This 
behavior is illustrated in Fig. [| As we have noted, this is important because the a' 
state couples to Standard- Model fields with an unsuppressed coupling 1//pq rather 
than with a volume-suppressed coupling 1//pq. 

Note that as y — > oo (reproducing the four-dimensional limit), the time needed 
for the probability to drop to zero increases without bound. Thus, in the four- 
dimensional limit y — > oo, the probability begins and remains at 1, as expected. 
However, for all finite values of y, the axion probability P a r-> a >(t) drops to zero in 
finite time, and remains there (as shown in Fig. |^). Indeed, for n max ^> 1, we see 
that the time To for the axion probability to drop to a predetermined fraction of its 
initial value scales as[] 

ro ~ -f- . (3.15) 

77 

' "max 

* Note that the derivations of these scaling results can often be quite subtle. In order to derive 
(3.15), we first observe that tq ~ to/TV 2 , where N is the normalization factor in (2.14) and where 



to is a time scale that is independent of n max . This scaling behavior holds independently of y 
when n max 3> 1, as can be verified numerically. Given this, it is straightforward to investigate the 
behavior of to as a function of y, leading to the result to ~ y 2 - We thus obtain To ~ {n max y)~ 2 . 
The final step is to realize that there is a hidden y-dependence buried in the meaning of the cutoff 
n max , and that in order to compare cutoffs for different values of y, we must choose a uniform 
y-independent convention for the Kaluza-Klcin truncation. Specifically, for each value of y, we must 
choose an appropriate y-dependent normalization of the cutoff n max such that the rescaled cutoff 



n max nas a fixed, y-independent net effect in an eigenvalue sum such as (3.9). It turns out that 
such a renormalization compels us to choose n' max ~ y 2 n max . We therefore find r ~ y 2 /(n' max ) 2 . 
To provide some explicit numbers, let us define tq as the time needed for the axion probability 
to fall to 10% of its initial value, and let us likewise define ?7 cr it to be the minimum number of 



Kaluza-Klein eigenvalues that must be included in the sum in (3.9) order to produce an initial axion 
probability of 0.99. (In this connection, note that it is only in the formal n max — > oo limit that 
the initial probability truly approaches 1.) We then find that to ~ 9.665/y 2 and n cr ; t » 981/y 2 . 
This latter relation enables us to normalize our values of n max in relation to 77 C rit by defining 
n max = "max/^crit = 2/ 2 «max/981, leading to the final result 

* 10-5 ufv- ( 3 - 14 ) 

I "max/ 

The value of the rescaled cutoff ri max is then arbitrary, and may be chosen according to considerations 
beyond those of our effective field-theory approach (such as truncating according to the underlying 
string scale). 
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Thus, as y — ► oo, we see that r — > oo, while for finite values of y this "decoherence 
time" is extremely short. This results in an essentially immediate suppression of the 
axion probability. Moreover, as we shall verify later in this section, this decoherence 
time is substantially smaller than the lifetime of the heaviest Kaluza-Klein mode 
contained within a'. This decoherence mechanism therefore renders the axion super- 
position a' virtually invisible with respect to subsequent axion interactions involving 



a'. 



It is straightforward to understand this suppression at an intuitive level. Unlike 
the case with the simple axion zero- mode do, in the present case our initial axion 
state is the infinite linear superposition a' given in ( |2.13| ). Let us therefore consider 
the behavior of the individual terms in the probability sum (|3.9| ) as n max gets large. 
For this purpose we may drop all factors of two and focus only on the behavior of 
Pa'^a' a s a function of n. max . For fixed y and large eigenvalues A, we have X 4 A\ « 2. 
Therefore, as we introduce increasingly heavy eigenvalues into the probability sum 
(|3.9|) , we find the effective behavior 



max 



E(1) + 2E(1)cos (M!f 

A A'<A V j 



(3.16) 



where each A-sum contains n max terms. For t = 0, this result factorizes to take the 
simple leading form 

^w(O) ~ J- ^(1 + 1) ~ -L . . (3.17) 

"max A, A' "max 

In other words, the presence of two independent A-sums provides an effective factor 
of nL a ^. which cancels the factor of n~5L that resulted from the normalization of a'. 
This initial coherent contribution of the independent two A-sums enables the initial 
probability to start at 1 regardless of the size of the volume of the compactified space. 
Indeed, the initial state a' is a highly coherent state. However, at later times t > 0, 
the cosine terms in ( |3.16D no longer add coherently to the sum, and their destructive 
interference effectively causes the sum to scale only as n max , corresponding to a single 
diagonal A-sum: 

Pa'-+a>(t) 1- ^(1) 1- ■ n max — . (3.18) 

"max A "max "max 

This phenomenon is completely analogous to the fact that a random walk traverses 
only the square root of the distance traversed by a coherent, directed walk. It is 
for this reason that the net axion preservation probability P a i^ a i(t) is so strikingly 
suppressed. Essentially, the initial axion state a' has "decohered" as a result of the 
incoherent Kaluza-Klein oscillations induced by the non-diagonal axion mass matrix. 
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We thus conclude that in higher dimensions, all axion-mediated processes which 
rely on the production and subsequent detection of the a' mode are strongly sup- 
pressed by a rapid "decoherence" which renders them virtually "invisible" . This de- 
coherence arises as a result of the destructive interference of the infinite-component 
laboratory axion oscillations. We see, then, that this provides an entirely new higher- 
dimensional mechanism which can contribute to the invisibility of a'-mediated pro- 
cesses. 

It is important to stress that this decoherence mechanism is relevant only for 
those processes which are sensitive to the time-evolution of the axion. As we have 
indicated, this includes all axion-mediated process {e.g., axion-exchange processes), as 
well as processes in which the axion is directly detected in the laboratory. Moreover, 
measurements of the axion flux from the Sun or from supernovae also fall into this 
category. However, this does not include processes which are insensitive to axion 
time-evolution. These include, for example, axion-production processes in which the 
axion appears only as missing energy. Nevertheless, invisibility can be achieved for 
such processes by adjusting the value of /pq via the mechanisms discussed in Sects. 2 
and 4. 

Finally, we remark that our higher-dimensional decoherence mechanism is com- 
pletely general, and applies not only to axions, but also to any bulk field <fi whose 
"shadow" interaction with the Standard-Model brane involves a coupling between 
Standard-Model fields and a "brane shadow" Kaluza-Klein superposition 

<P(y = o) = E c ^fc ( 3 - 19 ) 

k 

with non-zero coefficients c^. Just as in the axion case considered above, the non- 
trivial time-evolution of the Kaluza-Klein modes will cause the initial superposition 
(p(y = 0) to "decohere" extremely rapidly. This will necessarily produce a severe 
damping of any Standard-Model process on the brane that involves couplings to the 
coherent state <j)(y = 0). For example, this might therefore provide a partial solution 
to the notorious dilaton problem in string theory: it may simply be that the string- 
theoretic dilaton is "invisible", in much the same way as the axion is invisible. Similar 
considerations may also apply to Kaluza-Klein gravitons as well as other bulk moduli 
fields. 

3.3 Laboratory oscillations inducing a k — > a' 

Finally, let us consider a third relevant oscillation probability. As we discussed 
above, fields on the Standard-Model brane can couple only to the linear combination 
a'. This is why the probability P a '^ a '{t) is the appropriate probability for processes 
involving both Standard-Model production and detection of axions. However, for 
axions that are produced through mechanisms involving bulk fields (which are non- 
localized) rather than brane fields (which are localized), it is possible to envisage 



20 



situations in which a single Kaluza-Klein mode at {e.g., the zero-mode ao) is pro- 
duced. However, detection of such an axion mode on our Standard-Model brane 
continues to involve couplings to a', and therefore in such cases a relevant probability 
for detection is given by 



1 

N 
1 

N 



J2 A 4 A\ e~ Txt + 

A 

2 £ \ 2 (\') 2 AlAl,e-^ +r ^ 2 cos ~^ )2]t N 



A'<A 



(3.20) 



However, just as in the previous case, this probability vanishes in the limit n max — > oo, 
even if we again set T A = 0. Specifically, the probability to produce the specific 
coherent state a' on the Standard-Model brane vanishes as n max — > oo. Moreover, 
the same is true for all probabilities P afc _>„/(£) for all k and for all y. Of course, 
we cannot interpret this result as a decoherence, since even the initial probability 
at t — vanishes as n max — > oo due to the negligible overlap between and a'. 
Nevertheless, this demonstrates that even if the axion is produced by bulk fields not 
restricted to the Standard-Model brane, the probability for its subsequent detection 
in the laboratory is vanishingly small. 

3.4 Axion detection and decay 

Finally, we shall conclude this section by discussing a number of additional effects 
that arise due to the existence of an infinite tower of excited Kaluza-Klein axion 
states. Our comments in this subsection will primarily be focused on the possibility 
of axion detection and in particular on the role of axion decays. 

First, as already mentioned above, we must distinguish between processes in which 
the time-evolution of the axion plays a role and processes which are insensitive to 
the time-evolution of the axion. For example, the latter include axion-pro duct ion 
processes in which the axion is emitted into the bulk and therefore is manifested on 
the Standard-Model brane only as missing energy. A simple example of this is the 
axion-emission process 

FF -> a' -> bulk (3.21) 

where the axion, once produced, flies into the bulk. As we have seen in Sect. 2, this 
process scales as 1/ /pq rather than 1//pq. This already leads to the severe constraint 
/pq ~ O(TeV), and is completely analogous to the possibilities for detecting graviton 
emission in upcoming TeV-scale collider experiments. We shall discuss one such 
missing-energy signature in more detail in Sect. 5. 

By contrast, processes which involve an actual detection of the axion on the 
Standard-Model brane are necessarily time- dependent because they involve a non- 
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zero time interval between axion production and axion detection. However, even 
within this category of time-dependent processes, there are further subdivisions that 
can be made. One important distinction is the length of the time interval between 
axion production and axion detection. Certain axion-mediated processes, such as 
those taking place entirely within accelerator experiments, take place on time scales 
that are very short compared to the axion lifetime. By contrast, others (such as those 
involving fluxes of axions which are produced in the Sun or in supernovae and which 
are subsequently detected on Earth) involve much longer time scales. 

This distinction between "long" and "short" time scales therefore depends on 
the axion lifetime. Ordinarily, in four dimensions, the axion is relatively long-lived 
because it is so light and so weakly coupled to ordinary matter. Indeed, the dominant 
axion decay mode is to two photons, yielding the lifetime 

r 4D (a-77) = r « (— V f 4^ I « 10 48 seconds . (3.22) 

T 4D (a^77) \aj \m^ Q J 

In this expression, we are writing / P q ps 10 12 GeV to denote the usual four- 
dimensional Peccei-Quinn symmetry-breaking scale, and taking mpQ ~ 1CT 5 eV to 
denote the usual four-dimensional axion mass. (We choose these symbols in order to 
facilitate the comparison between the four-dimensional and higher-dimensional situ- 
ations.) We are also disregarding numerous model-dependent 0(1) coefficients which 
do not affect the overall scale of the result. Thus, in four dimensions, the axion is 
extraordinarily stable. 

In higher dimensions, this situation changes dramatically. Of course, the coherent 
axion mode a' is not a mass eigenstate, and thus, strictly speaking, it does not have 
a well-defined "lifetime" . Nevertheless, we can determine an effective lifetime for a' 
by estimating the shortest lifetime of any of the mass eigenstates a A of which it is 
comprised. This then yields the time scale over which the coherent state a' naturally 
decays as a result of its couplings to ordinary matter. In general, it is the heaviest 
Kaluza-Klein mass eigenstates which have the shortest lifetimes, with the decay mode 
into two photons continuing to be dominant. As a function of the cutoff n max , the 
"lifetime" of the coherent state a' can therefore be estimated to be 

r D>4 (a -> 77) = 




1(T 3 seconds . (3.23) 



In this expression we have taken -RMstrmg, implying A max w M str i ng . We have 

also chosen M str i n g = 1 TeV for simplicity. Thus, we see that the coherent state can be 
expected to decay to two photons much more rapidly than the usual four-dimensional 
axion. Of course, in this calculation we have taken the cutoff A max ps M string , which 
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represents the "worst-case" scenario. In any axion-related process of total energy 
E < M string , the production of Kaluza-Klein axion modes a\ with A > E will be 
kinematically disfavored. We would then take A max ~ E, which can increase the 
lifetime considerably. 

A priori, this significant reduction in the axion lifetime relative to the four- 
dimensional case means that our "decoherence" phenomenon for a' is irrelevant unless 
the decoherence time To is even smaller. However, it is straightforward to verify that 
this is indeed the case. Consulting ( |3.14j ), we see that 



Tq ~ y * ( P-) ( M^l) x 1(T 33 seconds (3.24) 
° y VI TeV/ VI TeVJ v ' 

where we have taken n' max « RM str i ng . For all phenomenologically interesting values 
of p and M s t r i n g, we thus conclude that r <C td>4- We see, therefore, that the 
coherent axion mode a' indeed decoheres sufficiently rapidly to justify the = 
approximation that was used in deriving it. [J] 

Given that the coherent state a' has an intrinsic lifetime r a i ~ 10~ 3 seconds, we 
can therefore use this as a benchmark for separating "long" and "short" processes. 
For laboratory oscillations over time scales shorter than this, it is legitimate to neglect 
the axion decay widths in calculations of these oscillations. By contrast, calculations 
of oscillations over longer time scales require the inclusion of the decay widths, and 
will therefore start to feel the effects of axion decays. 

It is therefore important to understand the effects of such axion decays, particu- 
larly as they relate to axion "invisibility". For this purpose, we may draw another 
distinction, this time between processes that take place entirely within a single de- 
tector (such as an axion-mediated process EE — > a' — > EE), and those which take 
place largely outside our detector (such as an axion beam travelling from the Sun 
to the Earth). For processes taking place entirely within a single detector, axion 
decay represents a breakdown of invisibility because we can in principle detect the 
emitted decay products. For example, even though we have found that P a /_» /(t) is 
suppressed for all times t exceeding the decoherence time To, the decohered state will 
nevertheless continue to propagate until the individual Kaluza-Klein axion modes 
that comprise this decohered state themselves decay. This will be discussed in more 
detail below. Given this observation, one might initially doubt the phenomenological 
importance of the a' decoherence phenomenon. However, the important point is that 
it is only the coherent state a' which couples to Standard-Model fields with the po- 
tentially dangerous unsuppressed coupling 1/ f P q, whereas individual Kaluza-Klein 
axion modes instead experience the safer suppressed coupling 1//pq. The deco- 
herence phenomenon therefore indicates that there is only a vanishingly small time 



'''In comparing ( |3.23 ) and (3.24), we must actually account for the relative rescaling of the 



cutoff n max — > n' mayl , as discussed above (3.14). This introduces an additional multiplicative factor 
(y 2 /981) 3 into ( |3.23 ), thereby shortening T£»4 by an additional factor « 10~ 9 for y w 0(1). 



However, we see from (3.24) that this still does not affect our main conclusion that tq <C to>4- 
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interval during which a process of the form FF — > a' — > FF can possibly occur with 
a dangerously large amplitude scaling as 1//pq. After this initial time interval, the 
a' state decoheres, and all subsequent interactions and decays will have amplitudes 
scaling as l//p Q . 

In sharp contrast are processes which largely take place outside our detectors. 
In such cases, we do not expect to be able to detect the decay products that result 
from axion decays at intermediate times. For example, in the case of an axion beam 
travelling from the Sun to the Earth, the photons emitted in flight through axion 
decays are presumably indistinguishable from the general background radiation. In 
such situations, therefore, axion decays lead not to a loss of invisibility, but rather 
to an enhancement of it. For example, let us consider the probability that the initial 
state a' will be found in a particular mass eigenstate a\ as function of time. Following 
the same procedure as above and using fl2.37 ), we find the probability 



P a ^a x (t) = IaM^"™. (3.25) 

Note that this is an exact result valid for all times, with axion decays producing an 
exponential suppression for the probability. Moreover, as we shall now demonstrate, 
this fact can be used to provide a direct experimental test of the higher-dimensional 
nature of the axion. Using ( |3.25| ), we may define 

P tot (t) £P^a A = ^E^^^ ( (3.26) 

A A 

as a "collective" amplitude that the axion mode a' survives as a function of time. This 
interpretation is justified because Ptot(^) is nothing but the time-dependent norm of 
the original a' superposition: 

Ptot(t) = (a'(t)\a'(t)) . (3.27) 

Given this norm, the collective decay width (i.e., the instantaneous decay probability 
per unit time) is given as 

1 dPtotit) J2\ Ta A 4 A\ e~ Fxt 

in complete analogy with the formalism for radioactive decays. However, because of 
the presence of an infinite tower of Kaluza-Klein states, it is immediately apparent 
that (T) is itself a function of time! For example, at very early times we have 
(r) ~ r Amax ~ M, 3 tring //| Q , where A max « M s t ring is the heaviest mass eigenvalue 
included in the linear superposition. By contrast, at extremely late times, we have 
(T) — > T\ where Ao is the lightest mass eigenvalue. Since (T) is the instantaneous 
decay probability per unit time (which can be measured in our Earth-bound detector), 
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a time-variation in (r) from different axion sources at different distances would serve 
as a direct experimental test of the higher-dimensional nature of the axion. 

Finally, let us briefly comment on the possibility of measuring the laboratory 
axion oscillations discussed previously. Although we have shown that probabilities 
such as Pa k ->a e (t) experience sinusoidal oscillations as functions of time, at a physical 
level these calculations presuppose that we are capable of experimentally detecting 
individual Kaluza-Klein modes a^. In other words, these calculations presuppose 
that we are able to distinguish the final axion states according to their Kaluza-Klein 
quantum numbers. In situations where such detections can be made, our previous 
results continue to apply. However, in the most straightforward scenarios, all Kaluza- 
Klein states will have identical decay modes {e.g., into two photons), leaving us with 
no experimental "handle" through which to detect the presence of an oscillation. 
Thus, the prospects for detecting laboratory axion oscillations depend crucially on 
the ability to perform laboratory axion measurements which are sensitive to particular 
Kaluza-Klein quantum numbers. Of course, this is completely analogous to the case 
of neutrino oscillations, where the existence of neutrino decay modes that distinguish 
between different neutrino quantum numbers (such as SU (2) gauge charge or flavor) 
permit the detection of neutrino oscillations. 

4 Cosmological relic axion oscillations 

In this section, we shall discuss a third higher-dimensional effect that can con- 
tribute to the "invisibility" of axions: the rate at which the energy trapped in cos- 
mological relic axion oscillations is dissipated. As we shall see, the presence of an 
infinite tower of Kaluza-Klein axion states can, under certain circumstances, actually 
enhance the rate at which this oscillation energy is dissipated. In such cases, the effec- 
tive Peccei-Quinn scale /pq can therefore be raised beyond its usual relic-oscillation 
bounds, leading to an axion whose couplings to matter are even more suppressed 
than in four dimensions. This can therefore provide another higher-dimensional way 
of achieving an "invisible" axion. 

Let us first recall the situation in four dimensions. One of the most important 
constraints on the scale of Peccei-Quinn symmetry breaking, and hence on the mass 
of the axion and the strength of its couplings to ordinary matter, comes from cos- 
mological relic axion oscillations. Unlike the laboratory oscillations discussed in the 
previous section, these cosmological relic oscillations arise due to the fact that as the 
universe cools and passes through the QCD phase transition at T m Aqcd, instanton 
effects suddenly establish a non-zero axion potential where none previously existed. 
In the usual four-dimensional situation, the axion can therefore find itself displaced 
relative to the newly-established minimum of the potential, and begin to oscillate 
around it according to the differential equation 

d (x _ _ r , . dct '~i ^ / . \ 

— + 3H{t) — + m 2 a a = , t> t QCD . (4.1) 
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Here H(t) = 1 / (2t) is the Hubble constant, where we are assuming for simplicity that 
this time-evolution takes place entirely in a radiation-dominated universe. We are also 
assuming that the axion mass m a is independent of time for t > £qcd and vanishes 
for times t < £qcd- We are also neglecting the decay width of the axion. In general, 
the initial amplitude of this oscillation at £qcd is set by the initial random angular 
displacement of the axion field. This scales as /pq, where /pq is the axion decay 
constant. Although these oscillations are ultimately damped due the cosmo logical 
Hubble expansion term, a relic of these axion oscillations should still exist today. 
Imposing the requirement that the energy stored in the relic oscillation today be 
less than the critical energy density (so as not to overdose the universe) then sets 
an upper bound /pq < 10 12 GeV which is consistent with bounds obtained through 
other means. This is therefore a bound on the "invisibility" of the axion. 

At first glance, it might seem that this issue should no longer play a role in 
our higher-dimensional scenario in which the fundamental scale /pq of Peccei-Quinn 
symmetry breaking is substantially lowered (perhaps even to the TeV-range) as a 
result of the volume factor in Q2.11| ). Indeed, it is the (low) fundamental mass scale 
/pq rather than the (high) effective mass scale /pq which sets the size of the VEV of 
the axion field a in the five-dimensional Lagrangian Q2.6|) , and which similarly sets the 
overall scale of the initial random angular displacement of the axion field. However, 
in passing from (|2.6| ) to the four- dimensional Lagrangian ( |2.9|) , there is an implicit 
volume-dependent rescaling of the axion field, as discussed below ( |2.1(J| ). Thus, even 



though initial displacement of the unrescaled axion field scales with /pq, the initial 
displacement of the rescaled axion field scales with /pq. The danger of an excessive 
oscillation energy density today is therefore just as relevant for our higher-dimensional 
scenario as it is for the usual four-dimensional case. 

In fact, this danger may actually be greater in the higher-dimensional case. This 
is because the Kaluza-Klein reduction yields not only the usual zero- mode axion do, 
but also an infinite tower of excited Kaluza-Klein axions at (k > 0). The above 
differential equation for the axion oscillations then generalizes to 

^ + 3if(t)^ + Mlta e = 0, t > t QCB (4.2) 

where A4 2 is the non-diagonal mass matrix given in ( £271 ). (Note that we shall 
continue to neglect axion decay widths in these equations; this assumption will be 
justified at the end of this section.) Of course, due to their Kaluza-Klein masses, these 
excited Kaluza-Klein axions feel a non-zero potential even prior to the QCD phase 
transition (i.e., even prior to the "turn-on" of ttt-pq), and it is therefore reasonable 
to assume that they are already sitting at their minima at the time of the QCD 
phase transition. However, due to the non-diagonal mass mixing matrix in (|4.2p , the 
initial displacement of only the zero-mode axion ao is sufficient to trigger the excited 
Kaluza-Klein modes into oscillation. This situation is illustrated in Fig. [|. 

Given this observation, there are a priori three possible effects that these excited 
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Figure 4: A plot of the coupled Kaluza-Klein cosmological relic axion oscillations. For this 
plot, we have normalized the initial displacement of the axion zero- mode to 1, and taken 
y = 0.5, m-pQ = 2, and £qcd = 1 (in dimensionless units). We have also considered the 
effects of only the first three excited Kaluza-Klein modes. Although the excited Kaluza- 
Klein modes have vanishing initial displacements, they are triggered into oscillation as a 
result of the initial displacement of the zero-mode. This in turn changes the subsequent 
time-evolution of the zero-mode. By contrast, the superimposed dashed line shows the 
behavior of the usual four-dimensional axion zero-mode in the y — > oo limit (when no 
Kaluza-Klein modes are present). 

Kaluza-Klein states can have on the system. First, it is possible that these excited 
Kaluza-Klein states will "capture" oscillation energy from the zero-mode oscillation, 
and essentially store it. Thus, in this case, the total energy density of the system 
would dissipate more slowly, thereby resulting in a greater relic energy density today. 
This would seriously strengthen the five-dimensional axion bounds relative to the 
usual four-dimensional bounds, and provide the most serious threat to the viability 
of these higher- dimensional scenarios. The second possibility is that although the 
excited Kaluza-Klein states steal energy from the zero-mode oscillation, they may 
be able to dissipate it more effectively. This would then lead to an enhanced relic 
energy loss rate, implying a weakened bound on the higher- dimensional scenarios. 
Finally, the third possibility is that the two effects cancel exactly, with the excited 
Kaluza-Klein states oscillating in such a way that even though they capture some 
energy from the zero-mode, they also alter the time-development of the zero-mode 
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in an exactly compensatory manner. Thus the traditional four-dimensional bounds 
would remain unchanged. 

In order to determine which of these possibilities is realized, we must solve the 
coupled differential equations fl4.2f) . In so doing, we shall make the following assump- 
tions. First, we shall assume that prior to the QCD phase transition at £qcd, the 
axion zero-mode experiences no potential and therefore has an initial displacement 
set to 1 (in units of fpq). We shall likewise assume that da /dt — at t = £qcd- By 
contrast, for the excited Kaluza-Klein modes (k > 0), we shall begin with the ini- 
tial conditions = da^jdt = at £qcd- These conditions reflect the fact that these 
excited modes have non-zero Kaluza-Klein masses even prior to the "turn-on" of mpQ 
at £qc;d, and hence should have essentially settled into their minima prior to the QCD 
phase transition. Although it is possible for the lighter Kaluza-Klein axion modes 
to have small displacements as well, for simplicity we shall ignore this possibility in 
what follows. Note that this assumption also enables us to start with identical ini- 
tial relic oscillation energy densities in both the four- and five-dimensional situations, 
and thereby enables us to make a direct comparison of the effects of the Kaluza-Klein 
modes on the cosmological time-evolution of the system. Finally, in substituting the 
matrix (|2.27|) into (|4.2j ), we shall assume that mpq(t) is given exactly by the step 
function 

m PQ (t) = mp Q e(t-t Q cD) • (4.3) 

The constancy of mpq(t) for t > £qcd is often referred to as the "adiabatic" ap- 
proximation. We caution, however, that for large values of mpQ this step-function 
approximation can differ quite substantially from the results of a more careful analysis 
of the time/temperature-dependence of the axion mass due to instanton effects [0. 

Given these assumptions, it turns out to be possible to solve the differential 
equations ( |4.2| ) analytically for all times t > £qcd and for an arbitrary number of 
Kaluza-Klein modes. We therefore do not need to make any further approximations 
(such as the traditional separation into so-called "overdamped" , "underdamped" , and 
"critically damped" oscillation phases). The first step in our analytical solution is 
to decouple the differential equations ( |4.2| ) by passing to the mass-eigenstate basis 
d\ defined in (|2.31|) . In order to work with dimensionless quantities, we shall define 



a\ = d\/ fpQ. Each of our uncoupled differential equations then takes the form 

^ + ff + ^-°- <*■*> 

For simplicity we can recast this equation into the form 

d 2 d x 3 dd x „ 

~*pr + Yr^ + ax = (4 ' 5) 

where we have defined r = At. The most general solution to this equation is then 
given by 

a(r) = r- l ' A [ C+ J 1/4 (r) + c_J_ 1/4 (r)l (4.6) 
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where J v (t) are the Bessel functions of first kind. We therefore wish to solve for the 
unknown constant coefficients c + and c_. To do this, we impose our initial conditions. 
In the original Kaluza-Klein basis, these conditions are given by afc(ro) = /pq<5/co 
and dak(To)/dr = where To = tqcd is the initial time at which we begin the 
time-evolution of our axion fields (representing the "turn-on" time for ttipq). In the 
dimensionless mass-eigenstate basis, these initial conditions therefore take the form: 



a A (r ) = A x , 



dr 



(4.7) 



T=T 



where A\ is defined in ( |2.33| ). The assumptions underlying these initial conditions 
were discussed above. Thus, in the mass-eigenstate basis, we see that each of the 
Kaluza-Klein modes begins with an initial displacement. Solving for c± is then 
straightforward. The first initial condition (the displacement condition) trivially 
gives the constraint 

c + J 1/4 + c_J_ 1/4 = A x tI /4 (4.8) 

where for notational convenience any Bessel function written without an argument 
is understood to be evaluated at r . Using the Bessel-function identity 

j^ ±V Ur)\ = ±r^J UTl (r) } (4.9) 
we see that the initial velocity constraint takes the form 

c + J 5/4 (r) + c_J_ 5/4 (r)l , (4.10) 



dd\ _— 1/4 



T 



dr 

implying c+Js/4 = c_ J_5/4. Together with (|4.8|), this leads to the solutions 



c± = -^A x r* /4 J T5/4 (4.11) 
where we have used the further identity 

J1/4J-5/4 + J-1/4J5/4 = • (4-12) 

TTTo 

Substituting c± from ( |4.11| ) into the solution given in (fOJ), we thus obtain our 
final closed-form solution to the axion differential equation: 

2a(t) = —A^V^j^r) (4-13) 

where we have defined 

j(r ;r) = J- 5 /4(t )J 1/4 (t) + J 5/4 (t )J_ 1/4 (t) . (4.14) 
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This implies that the first time-derivative is given by 

dr ~ V2 X ° } 
where we have likewise defined 

/(to;t) = J-s/4,(t ) J 5 /4(t) - J5/i(r ) J_ 5 / 4 (r) 



(4.15) 



(4.16) 



Note that j'(to; r) — ► as r — ► To, as expected, since the initial velocities vanish for 
each of the axion modes. 

Now, the energy contribution from a single mode d\ is given by 



Px(r) = y 



~2 , ( da * 



(4.17) 



where A is the dimensionless eigenvalue defined in ( |2.29| ) and where p = pj (^pq/pq) 
is a dimensionless energy density. Substituting in the above results, we therefore find 



7T 



(4.18) 



While this expression is exact for all times r, it is also useful to have an ap- 
proximation valid for extremely late times satisfying r ^> 1 (such as the present 
cosmological time). Using the Bessel-function asymptotic expansion 



JJt) PS \ COS T - 

TIT V 2 4 



as t — > 00 



(4.19) 



it is straightforward to show that 



2 

./(to; r) 2 + j'(r ; r) 2 ps — {[J 5 /4(to)] 2 + [J- 5 /4(r )] 2 + V2 J 5/a (t )J_ 5/4 (t )} . 

(4.20) 

Given the exact result Q4.18 ), we can convert from r back to our original time 
variable t to obtain the final closed-form solution 



7T 



p x (t) = y A\ A 4 t 5/2 r l/2 [j(At ; At) 2 + j'(At ; At) 2 



(4.21) 



where we have defined the dimensionless time t = m^qt. Thus, adding together the 
contributions from all the mass-eigenstate modes, we obtain the final energy density: 



~pif) = EpaW = \ t 5/2 t- 1/2 E^A 4 [j(Xto;\t) 2 +fCXt ;\tf 

A A 



(4.22) 
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Note that this is the exact result for the relic oscillation energy density as a function 
of time. 

It is easy to verify that pit) has the correct limit as t — > to- Indeed, as t — * to, we 
find using the identity Q4.12Q that pit) — > p , where 



Po^£AM*. (4.23) 

This is indeed the correct value of the initial energy in the mass-eigenstate basis, 
since each mass eigenmode starts with zero velocity and with initial displacement 
A\. However, in our original Kaluza-Klein basis, our initial conditions at to consist 
of having only the zero- mode displaced by 1. This implies that p — 1/2. Comparing 
this result with (|4.23|) then yields the identity quoted in ( |2.35| ). In particular, this 
identity holds for all y, as can be verified directly by substituting the values of A and 
A\ and evaluating the eigenvalue sum. 

For late times At ^> 1, the exact result ( |4.22|) for the energy density simplifies to 
take the form 

p(t) = ^X(to)t* /2 t-^ (4.24) 
where the time-independent coefficient X{to) is given by 
X(t ) = £^A 3 {[J 5 /4(At )] 2 + [^5/4(^ . (4.25) 

A 

This demonstrates that asymptotically, the total energy density falls as p(t) ~ t~ 3 / 2 
regardless of the presence of the excited axion Kaluza-Klein modes. Moreover, in the 
"double-asymptotic" cases in which we also have Ato 3> 1 for all A, we can further 
approximate 

X(t ) « -iE^i = "i, (4-26) 

Tito A 7Tt 



where we have used the identity ( |2.35| ) in the last equality. We then find 



-3/2 

pit) « ~ ( f ] , (4.27) 




which is consistent with the initial energy density p(to) = 1/2. However, for practical 
purposes, we shall focus on the expressions ( [4.24| ) and ( |4.25| ) in which t is taken 
large but to is kept arbitrary. This is because Ato can occasionally be relatively small 
(particularly for the lightest eigenvalues) even if to itself is relatively large. Finally, 
in the four-dimensional y — > oo limit, we find that A^ — > ky for k > 1, while Ao — > 1. 
Thus, in this limit, the excited Kaluza-Klein modes decouple and we obtain the 
four-dimensional result 

X 4D (to) = [</ 5 /4(t )] 2 + [J-5/4 (to)? + V2J 5/4 (to)J- 5 /4(to) ■ (4.28) 
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As expected, this approaches the value given in ( |4.26| ) for t ^> 1. 



We see from ( |4.24|) that although the rate of energy loss remains fixed at asymp- 
totic times, the presence of the Kaluza-Klein states can nevertheless change the overall 
value of the energy. This is because the oscillating excited Kaluza-Klein states can 
change the rate of energy loss at intermediate times, possibly leading to an enhanced 
or diminished energy at late times. The cumulative effect of these Kaluza-Klein states 
at late times is encoded within the expression X(t ). Therefore, in order to under- 
stand the effect of the Kaluza-Klein states on the energy dissipation rate, we need to 
understand the behavior of X(t ) as a function of the radius variable y = (mpqi?) -1 
for a fixed initial time to. 

The results are rather surprising. Of course, for to ^> 1, we are in the "double- 
asymptotic" regime (t ,t) ^> 1 for which ( |4.26| ) and fl4.27p are expected to apply. 



We therefore find that in such cases the presence of the Kaluza-Klein modes does 
not alter the energy density relative to the energy density that would have been 
obtained in four dimensions. In other words, even though we have an infinite set of 
Kaluza-Klein axion modes which are induced into oscillation as a result of the initial 
displacement of the axion zero-mode, these oscillations nevertheless change the time- 
development of the zero-mode in a compensatory manner so that the total oscillation 
energy density as a function of time is exactly preserved. This indicates that in such 
situations, the higher-dimensional axion scenarios are no less viable than the usual 
four-dimensional scenarios. 

Even more surprising, however, is the situation that arises for smaller to- I n 
such cases, we cannot use the "double-asymptotic" expression ( }4.26| ), and we must 



resort to the full expression in ( |4 . 2 5| ) . We then find that X(t ) is smaller than the 
four-dimensional value given in (|4.28| ), which implies that 

p(t) < pm{t) - (4.29) 

In other words, even though the excited Kaluza-Klein states are triggered into oscilla- 
tion by the initial displacement of the zero-mode, these Kaluza-Klein states dissipate 
the oscillation energy more efficiently and result in a net decrease in the oscillation 
energy as a function of time relative to the four-dimensional case. This situation 
is shown in Fig. || where we plot pit) / p^it) = X(to)/X 4 D(to) as a function of 
y = (mpQ_R) _1 for three different values of to- Note that in the y — > 0, to — > limit, 
the ratio pit)/ p^it) decreases as ^Jy. 

We thus conclude that in such situations, the presence of coupled relic Kaluza- 
Klein axion oscillations can actually weaken the usual four- dimensional upper bounds 
on fpQ. This implies that it may be possible to consider higher values of /pq than 
are usually allowed in four dimensions, thereby further diminishing the axion cou- 
plings to matter and providing yet another higher-dimensional method of achieving 
an "invisible" axion. We hasten to point out, however, that the size of this effect 
depends crucially on the value of y as well as on the initial time to at which the axion 
potential is established. In general, for a given value of to, we can expect to see a 
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Figure 5: The energy-dissipation ratio factor p{t) / p^£,(t) = X(to)/Xm(to) as a function of 
y = (mpQi?) -1 , assuming (a) to = 10; (b) to = 1; and (c) to = 0.1. All cases reduce to the 
the usual four-dimensional result as y — > oo, with the Kaluza-Klein axion states decoupling. 
However, as the size of the extra dimension grows and y decreases from infinity, we see that 
the net effect of the Kaluza-Klein states is to dissipate the relic oscillation energy density 
more rapidly, leading to smaller relic oscillation energy densities at final times. The size 
of this effect depends on to, with smaller values of to corresponding to sizable effects at 
large values of y, while for larger values of to this effect is delayed until correspondingly 
smaller values of y. In general, this effect becomes substantial for yto ^$ 0(1). Note that 
all curves tend to zero in the y — > limit, implying an infinitely rapid dissipation of the 
relic oscillation energy density in the full five-dimensional limit. 

sizable deviation from the four-dimensional asymptotic result only when Ato 0(1) 
for the lightest eigenvalues. This implies yto < 0(1). Taking t = 10~ 5 seconds 
(corresponding to the QCD phase transition) and a reference value mpQ = 10~ 4 eV, 
we find to = m PQ^o ~ 10 6 . (By contrast, the current cosmo logical time is given by 
t pd 10 29 , corresponding to t = 10 18 seconds.) This implies that in practice, we should 
not expect to see a sizable decrease in the cosmological relic oscillation energy density 
unless y < (9(10~ 6 ). However, this assumes that a particular value of m P q, which in 
turn assumes a particular fixed value of /pq via ( |2.25| ). Changing /pq can therefore 
change this result substantially. We shall discuss this possibility in Sect. 5. 

Finally, let us briefly discuss the issue of axion lifetimes as they relate to these 
cosmological relic oscillations. Ordinarily, relic axion oscillations are important pre- 
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cisely because the usual four-dimensional axion is so long-lived. However, in higher 
dimensions, we have seen in Sect. 3 that the heavy Kaluza-Klein axion modes be- 
come more and more unstable, and hence cannot be expected to survive over the long 
cosmological time scales we have been assuming. It is therefore natural to wonder 
whether this places our conclusions about enhanced energy dissipation rates in jeop- 
ardy. In other words, it is not a priori obvious that we are justified in neglecting the 
excited Kaluza-Klein decay widths in (|4.2| ). However, we have seen that relic axion 
oscillations in higher dimensions are primarily sensitive not to the heaviest axion 
modes, but to the lightest modes. For sufficiently large radii, these lightest modes 
have masses which are relatively close to the mass of the axion zero-mode, and we 
have already seen in Fig. [I] that this in turn is bounded from above by the usual four- 
dimensional axion mass mpQ. Thus, the lightest axion Kaluza-Klein modes continue 
to be extremely long-lived, and will therefore survive to induce the enhancement of 
the oscillation energy dissipation rate that we have observed. 



5 Some numbers, bounds, and constraints 

Let us now combine our different results from the preceding sections. Our goal 
will be to determine the extent to which a self-consistent picture of axion energy 
scales emerges from the previous results. 

We begin with the three fundamental equations given in ( |2.11|) and (|2.25|) . For 



simplicity, in Q2.25| ) we shall take g 2 = (?q UT = 1/2, and in ( |2.11| ) we shall take 
5 = 1. We can therefore combine these three equations in order to express y in terms 
of the fixed quantities Aqcd, M stT i ng (the fundamental underlying mass scale in the 
theory), and the effective axion decay constant /pq (our measure of "invisibility"). 
For simplicity we shall also take /pq = M str i ng , since we want to have only one 
fundamental mass scale in the problem. This then yields the result 

£ /pqAqcd 

We know that M string cannot be lower than « (9(1 TeV), and Aq C d is absolutely 
fixed (~ 250 MeV). We therefore find that in general, y is bounded from below: 

2.5 x 10 12 GeV , . 

V Z Z/min ~ — ~ • (5.2) 

£ JPQ 

Thus, assuming £ « 0(1), we see that /pq ~ 10 12 GeV is consistent with having 
y ~ 0(1). Remarkably, this is precisely the region where we expect to find the axion 
mass becoming independent of /pq, as shown in Fig. [l|. 

While this shows the self-consistency of the situation depicted in Fig. [I], a natural 
question arises as to whether it is possible to tolerate larger values of M string . Indeed, 
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slightly larger values of M str i ng [e.g., in the O(10 TeV) range] may be preferred on 
the basis of detailed comparisons with experimental data. Ordinarily, it might seem 
to be impossible to increase the value of M string any further, because we see from 
(|5.1|) that increasing M string requires increasing / P q in order to maintain the 0(1) 
values of y (as preferred on the basis of Fig. [I]), and increasing /pq generally runs into 
difficulties with cosmological relic oscillation energy densities overdosing the universe. 
However, we have seen in Sect. 4 that the Kaluza-Klein axion modes may be capable 
of dissipating this excess energy density more rapidly so as to evade these bounds. 
The question then arises: to what extent can we increase /pq, thereby making the 
axion increasingly "invisible" , without disturbing the relic energy density bounds? 

Note that increasing /pq has a number of effects. First, as /pq increases, we 
find from (|5.2|) that y m ; n decreases. This means that y can be chosen even smaller. 
However, we also find from ( |2.25| ) that mpQ decreases, which in turn implies that 
to = mpQ^o (the dimensionless time of the QCD phase transition) also decreases. 
Defining y cr - lt to be the critical value of y at which we start to observe a significant 
decrease in the relic oscillation energy density, we have already seen in Sect. 4 that 
Z/crit ~ £o • Thus, as t decreases, we see that y CTit increases, implying that it becomes 
easier to compensate for the effect of having increased /pq in the first place. Indeed, 
this suggests that there might be an alternative, self-consistent, significantly higher 
value of /pq than previously thought. 



In order to determine this self-consistent value of /pq, we first note from ( |2.25| ) 
that 

(2.5 x 1(T 3 GeV 2 K , v 

m PQ « i 11 . (5.3) 

/pq 



This in turn implies that 



(3.8 x 10 16 GeV) f 
mpQto ~ * — (5.4) 



where we have taken to ~ 10 5 seconds (corresponding to the QCD phase transition). 
We thus have 

/pQ 

^ crit ~ 3.8 x 10"> GeV ' (5 ' 5) 
which implies that we can obtain y crit « 0(1) simply by taking £ rs 0(1) and 

« 3.8 x 10 16 GeV w M GUT . (5.6) 



In other words, if we take /pq ~ Mgut ; then we have a self-consistent solution with 
Umm ~ 0(X), 2/crit ~ 0(1), and M string w 20 TeV. For such values, the axion mass 
is independent of /pq, and the Kaluza-Klein modes begin to induce a significant re- 
duction in the final relic oscillation energy density which can in principle compensate 
for the increase in /pq. Remarkably, this analysis suggests that /pq, the effective 
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Peccei-Quinn symmetry-breaking scale, may be related to Mgut, the effective four- 
dimensional GUT symmetry-breaking scale. 

Of course, there is still one constraint that we have not imposed in the above 
analysis: we have not restricted the size of the radius R of the extra spacetime di- 
mension. In principle, this is not a problem because the Standard-Model fields are 
restricted to a D-brane, and thus there are no bounds on the sizes of such transverse 
extra dimensions that can arise from Standard-Model processes. However, gravity 
is generally free to propagate into whatever extra dimensions exist, leading to the 
additional constraint R < O (millimeter). It is therefore important to understand 
how this additional constraint limits the above scenarios. We stress, however, that 
imposing this additional constraint relies on the assumption that gravity is indeed 
free to propagate in the extra dimensions. Recent ideas concerning "gravity localiza- 



tion" [H| have shown that this need not always be the case. 

If we do assume this to be the case, however, then the above scenarios are signif- 
icantly restricted. Requiring R^ 1 > 10~ 4 eV impliesf] that 

. 10- 13 GeV 
™pq Z • (5.7) 

y 

Moreover, using ( |5.3| ), we find that this implies that 

/ PQ < (2.5 x 10 10 GeV) £y . (5.8) 

This therefore sets an upper bound on /pq as a function of y. Combining this upper 
bound with ( |5.1|) , we thereby obtain the constraint 

Thus, the lower limit for y depends crucially on the model-dependent parameter £ 
and the value we choose for M string . Since £ reflects the PQ charges of the ordinary 
fermions, it is not unreasonable to assume that £ may be somewhat larger than 1. 
Taking M str i ng ~ 1 TeV therefore still enables us to have y ~ 0(1). However, despite 
this fact, we can combine ( |5.5|) with (|5.8| ) to show that 

y CTit < (6.6 x l(T 7 )y . (5.10) 



* Note that imposing this constraint is actually somewhat subtle, and depends on a choice of 
which variables to hold fixed. In the analysis in this section, we have been taking /pq and y as 
inputs, and treating M str i ng , R, and mpQ as derived quantities. Thus, with this convention, /pq 
is considered to be independent of R, while /pq (identified with Muring) is considered to be an 
independent quantity. Note that this procedure exactly mimics the situations [To| in which the 
GUT and Planck scales are lowered by extra spacetime dimensions: it is always the "measured", 
large, four-dimensional scale that is held fixed, while the reduced higher-dimensional scale is viewed 
as a function of R. We therefore continue this convention in the present case even though neither 
/pq nor /pq has been experimentally measured. 
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Thus the value of y in this case is always significantly larger than the critical value 
that would be required in order to reduce the relic oscillation energy density below 
its usual four- dimensional value. 

Of course, this does not disturb the self-consistency of this scenario. As a result of 
(|5.8|) , /pq may still be in the range that satisfies the usual four- dimensional bounds, 
and we have seen in Sect. 4 that the presence of the Kaluza-Klein axion states does not 
increase the final relic oscillation energy density relative to the four- dimensional case. 
Moreover, these axions continue to be virtually "invisible" against direct detection 
and/or subsequent interactions as a result of the decoherence effect discussed in 
Sect. 3. Thus this picture continues to be self-consistent, and continues to lead to an 
invisible axion. Furthermore, although we have restricted our analysis to the case of 
a single extra dimension, the corresponding constraints in higher dimensions may be 
significantly weaker. In any case, a much more detailed analysis is necessary in order 
to make the above numerical bounds more precise, and to determine whether further 
experimental constraints may be imposed. 

There are also other phenomenological constraints that may be imposed, partic- 
ularly constraints that are insensitive to axion time-evolution. Good examples of 
this would be axion missing-energy signatures, such as might arise from the decay 
K + — > ir + a'. In four dimensions, this process has a branching ratio which scales 
as fpQ, leading to a bound /pq > 10 4 GeV. In five dimensions, by contrast, this 
branching ratio scales as 

r>p / jf+ + n 1 V?, ( mK \ 1 

BK(K + —> 7r a ) ~ -j^— 2^ 1 j2 — (5.11) 

JPQ n=0 JPQ V J staring/ JPQ 

For ttik ~ 500 MeV, this then implies the constraint 



/pq £ Jv^l0 4 GeV. (5.12) 

y ^"string 

Taking /pq = M string then leads to the bound M string > 370 GeV, which is consistent 
with the idea of lowering the string scale to the TeV-range. Of course, we reiterate 
that a much more detailed analysis of these and other processes is necessary in order 
to sharpen these bounds and constraints. 

Finally, it may also happen that an axion process in higher dimensions exactly 
reproduces the four- dimensional result. As an example of this, let us consider an a!- 
exchange process at zero momentum transfer (e.g., FF —> a' — > FF). The amplitude 
for such a process is given by 

1 i CO 

A = — (a' a') ~ r m r n(a m a n ) (5.13) 

•'PQ /pQ m,n=0 

where (BA) denotes the propagator from state B to state A. Passing to the mass- 
eigenstate basis d\ via ( |2.31| ) and using the zero-momentum propagator (aydx) = 
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^a'a/A 2 , we find that this amplitude then takes the form 

1 1 00 1 

A = T2 2~ ^ \2 ^ r m r nU X mUxn = — • (5.14) 

iPQ m PQ A ^ m,n=0 JPQ m PQ 

Note that we have used (|2.36|) followed by (|2.35| ) in the final equality. However, we 
see that the final result is nothing but the amplitude that we would have obtained in 
four dimensions for an axion that couples with the usual four-dimensional coupling 
/pq and has the usual four- dimensional mass mpQ. Thus, in this case, we obtain no 
new bounds coming from such a'-mediated processes. This stands in stark contrast 
to the analogous case of graviton-mediated processes, from which one can generally 
derive stringent bounds on the radii of the extra spacetime dimensions. Of course, 
we stress that this result holds only for zero momentum transfer, and is likely to be 
different when sizable momenta are carried by the intermediate axion state. 



6 The Standard-Model dilaton 



The above considerations about placing a Standard- Model singlet field in the bulk 
are actually quite general, and transcend the specific example of the QCD axion. To 
illustrate this point, let us briefly consider the case of another conjectured particle, 
the Standard-Model dilaton. This particle is introduced into the Standard Model 
in order to restore the classical scale invariance broken by mass terms. The scale- 
invariant extension of the scalar sector of Standard Model is given by 



L = \{d^Df + (Z^) -V ((f),a) 



(6.i; 



where <fi is the Higgs field and where the Standard- Model dilaton field D, like the 
axion field a, is written in terms of a decay constant fo via a relation of the form 



D = f D exp(a/f D ) 



(6.2) 



By suitably choosing the parameters in the tree-level Higgs potential Vo(0, a), we can 
arrange (a) = 0. Consequently (D) = fo represents the mass scale at which dilata- 
tion invariance is spontaneously broken. After quantum corrections are included, the 



scalar potential can be written in the form |21 



V 



If 



■v 



where V = 



D 



y o (0,O) + y«(0,O)-A(0,O)ln 



D 



f. 



D 



(6.3) 



Here V^(<j>, 0) is the one-loop contribution to the effective potential and A(0, 0) is 
the divergence of the dilatation current. The presence of A(0, 0) in this expression 
breaks the scale invariance, and gives rise to a dilaton mass 



ni 



D 



4(A) 
fh 



(6.4) 
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where (V) = A/4 at the minimum of the potential (|6.3| ). Note that in the Standard 
Model, the heavy top-quark mass leads to (A) > 0. Thus, in order to change the 
sign of (A) , one requires additional heavy Higgs-boson contributions for the stability 
of the dilaton potential. 

Let us now consider what happens when the dilaton field propagates in a five- 
dimensional bulk and therefore has a Kaluza-Klein decomposition of the form 



y) = £ D n (s")cos^|) • 



(6.5) 



n=0 

The five-dimensional action for the dilaton then takes the form 



S 



d 4 xdyM s \\{d M Df-V{x)5{y) 



(6.6) 



where V is given in ( |6.3|) . Substituting (|6.5 ) into (|6.6|) and integrating over the fifth 
dimension then gives rise to the effective four-dimensional Lagrangian 



oo n 2 



cff 



n=0 



\R 2 



E T n D n 



D \n=0 



V Q ((/>, 0) + V^\<j>, 0) - A(0, 0) In I j- £ r n D n 

\JD n=0 J 



(6.7) 



where fo = ^j27rRM stI i n gfD and where we have canonically normalized the dilaton 
kinetic terms. The minimum of the effective dilaton potential therefore occurs at 



(Do) = fo 

(D n ) = for all n > 1 



(6. 



We can derive the dilaton mass matrix by considering the local curvature of the 
effective dilaton potential near its minimum. This is completely analogous to the 
axion case (|2.22| ), and gives rise to 



/ « a \2 _ n 2 4A 

( M D) nn > = -jpO n n> ~ —r. 



(6.9) 



Remarkably, this mass matrix has exactly the same structure as in the axion case 
(|2.24|) , and consequently the physical implications will be identical to those for the 
axion. Of course, this result is expected since the dilaton and axion are both Nambu- 
Goldstone bosons of a spontaneously broken symmetry, and consequently have similar 
couplings to the anomalous divergences of their respective currents. 
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7 Conclusions 



In this paper, we have studied some of the novel effects that arise when the 
QCD axion is free to propagate in the bulk of large extra spacetime dimensions. 
First, we found that under certain circumstances, the mass of the axion can become 
independent of the energy scale associated with the breaking of the Peccei-Quinn 
symmetry. Because this energy scale determines the couplings between the axion and 
ordinary matter, we now have the freedom to adjust the strength of these couplings 
without disturbing the axion mass. This can therefore provide a new mechanism for 
achieving an invisible axion. 

Second, we pointed out that in such higher- dimensional scenarios, the axion will 
typically experience laboratory axion oscillations which are completely analogous to 
neutrino oscillations. This is therefore a new and unexpected phenomenological fea- 
ture for axions that does not exist in the usual four- dimensional case. Moreover, we 
found that these laboratory oscillations can cause axions to "decohere" extremely 
rapidly. This is therefore a second higher-dimensional phenomenon that may con- 
tribute to an invisible axion. Moreover, this phenomenon arises for all non-zero radii. 

Third, we discussed the role that excited Kaluza-Klein axion states may play in 
axion-mediated processes and axion decays. This enabled us to propose several direct 
experimental tests of the proposed higher-dimensional nature of the axion. 

Finally, we found that under certain circumstances, the presence of these Kaluza- 
Klein axion modes can significantly accelerate the dissipation of the energy associated 
with cosmological relic axion oscillations. Moreover, even when these circumstances 
are not met, we found that the Kaluza-Klein states do not induce a violation of 
the usual four- dimensional bounds. This demonstrates that such higher-dimensional 
axion scenarios are no less viable than their four-dimensional counterparts, and indeed 
may even be preferred on the basis of their remarkable "invisibility" decoherence 
properties. 

Of course, there are many aspects of higher-dimensional axion phenomenology 
which we have not examined in this paper. These include the role that axions play in 
stellar evolution, the thermal production of axions, axionic string decay, and isocur- 
vature axion fluctuations. While some of these topics have been discussed in Ref. [jCJ . 
it will be interesting to further explore the role that extra spacetime dimensions can 
play in these areas. 

Although we have focused primarily on the case of QCD axions, we stress that 
much of our analysis is completely general and may apply for other bulk fields as 
well. This was explicitly illustrated in Sect. 6, where we considered the case of the 
Standard-Model dilaton. Similarly, we expect that our analysis will also apply to 
other bulk fields such as Kaluza-Klein gravitons, string-theoretic dilatons, and other 
bulk moduli. Indeed, the twin properties of laboratory oscillations and decoherence 
leading to "invisibility" are likely to play an important role in experimental searches 
for such particles, and likewise an analysis of the effects of their excited Kaluza- 
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Klein modes on cosmological evolution is likely to parallel the analysis in Sect. 4. 
Moreover, both of these effects are likely to play an important role in the all-important 
questions of dilaton and radion stabilization. We therefore leave these issues for future 
investigation. 
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